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Theoretical spectra and evolutionary models that span the giant planet brown dwarf con- 
tinuum have been computed based on the recent discovery of the brown dwarf, Gliese 229 B. 
A flux enhancement in the 4 5 micron window is a universal feature from Jovian planets to 
brown dwarfs. We confirm the existence of methane and water in Gl 229 B's spectrum and 
find its mass to be 30 to 55 Jovian masses. Although these calculations focus on Gliese 229 B, 
they are also meant to guide future searches for extra-solar giant planets and brown dwarfs. 
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Brown dwarfs inhabit a realm intermediate between the more massive stars and the less massive planets. 
Their thermal infrared emission is powered by the release of gravitational potential energy as regulated by 
their atmospheres. Long known only as theoretical constructs, the discovery of the first unimpeachable 
brown dwarf (1,2) allows a detailed study of a representative of this population of objects. Gliese 229 B, 
the recently-discovered companion to Gliese (Gl) 229 A, has an estimated luminosity of 6.4 ± 0.6 x 10~ 6 L Q 
(solar luminosity), an effective temperature, T e g, below 1200 K, and a clear signature of methane in its 
spectrum (3). Since there can be no stars cooler than 1700 K, with luminosities below 5 x 10~ 5 L Q , or 
with methane bands (4), Gl 229 B's status as one of the long-sought brown dwarfs is now beyond question. 
However, models of Gl 229 B's atmosphere and evolution are required to derive its physical properties and 
the previous lack of observations had inhibited the generation of theoretical spectra. To remedy this, we 
coupled model spectra and evolutionary calculations to estimate the object's T e g , L, surface gravity g, mass 
M, radius R, and age t, and to find useful spectral diagnostics. The recent discoveries of planets 51 Pegasi 
B, 70 Virginis B, 47 Ursae Majoris B, and Gl 411 B (5) have doubled the number of known Jovian planets. 
There is now an extraordinarily rich variety of low-temperature, low-mass (0.3 - 84 Mj (Jupiter mass)) 
planets and brown dwarfs. Our improved evolutionary models and spectra, here applied to Gl 229 B, are 
meant to facilitate the study and interpretation of these objects. 

To compute the atmospheric temperature profile for brown dwarfs in the relevant temperature range 
(600-1200 K), we adapt a model originally constructed to study the atmospheres of the Jovian planets and 
Titan (6). We assume a standard solar composition for the bulk of the atmosphere (7). Refractory elements 
(for example Fe, Ti, and silicates) condense deep in the atmosphere for T e s « 1000 K, and thus have 
negligible gas-phase abundance near the photosphere, as is also true in the atmosphere of Jupiter (T e s = 124 
K). For an atmosphere similar to that of Gl 229 B, chemical equilibrium calculations indicate that C, N, O, 
S, and P are found mainly in the form of methane (CH4), ammonia (NH3), water (H2O), hydrogen sulfide 
(H2S), and phosphinc (PH3), respectively. However, deep in the atmosphere, chemical equilibrium favors 
CO over CH4 and N 2 over NH 3 . Our model atmosphere incorporates opacities of these molecules, H 2 , and 
He (8) in their respective solar abundances and includes no other elements. 

To constrain the properties of Gl 229 B, we construct a grid of brown dwarf model atmospheres with T e s 

2 



ranging from 600 to 1200 K and 100 < g < 3200 m s~ 2 . For each case we compute a self-consistent radiative- 
convective equilibrium temperature profile and the emergent radiative flux (9) . Absorption of radiation from 
Gl 229 A is included in our model, but contributes negligibly to Gl 229 B's energy balance owing to the large 
orbital separation (> 44 AU) and faintness of Gl 229 A. 

Emergent spectra of brown dwarf atmosphere models compared to observed fluxes (Fig. 1) (1,10) show 
the influence of a minimum in the molecular opacities at wavelengths around 4 — 5 fim. As in the case 
for Jupiter, this minimum allows radiation to escape from deep, warm regions of the atmosphere. Clearly, 
this wavelength region is advantageous for future brown dwarf searches. By comparison, the widely-used K 
band at 2.2 fim is greatly suppressed by strong CH4 and H 2 — H 2 absorption features. Beyond 13 /jm, the 
decreasing flux falls slightly more rapidly than a Planck distribution with a brightness temperature near 600 
K. 

Our computed spectra (Fig. 1) are a good match with the data in the 1.2 — 1.8 ^m window regions, 
but deviate at 1 /im, in the window centered on 2.1 /zm, and in regions of low flux. Our best-fitting models 
reproduce the observed broad band fluxes (3) reasonably well. While many individual spectral features 
of CH4 and H 2 are reproduced, particularly near 1.7 and 2.0 /zm, the overall band shapes are not well 
accounted for in the 1 - 2.5 ztm region (Fig. la). We attribute these discrepancies to a poor knowledge of 
the CH 4 opacity and, to a lesser extent, the H 2 opacity. Although we have combined several sources of 
varying accuracy (8) to generate as complete a description of the CH 4 opacity as possible, methane line lists 
are based on laboratory measurements at room temperature and do not include lines from higher energy 
levels that would be populated at brown dwarf temperatures. Thus, the opacity of CH4 atTsj 1000 K is 
the most likely cause of the mismatches seen in the 1.6 - 1.8 /im band and at A > 2.1 ziin. 

Clouds may alter the atmospheric structure and spectrum of Gl 229 B, as they do in the atmospheres 
of planets of our solar system. Extrapolating from results for Jupiter (11) and using more recent chemical 
equilibrium calculations (12), we find that the following additional molecules are expected to condense 
between 10~ 3 and 10 bars: NH4H 2 P04, ZnS, K 2 S, Na 2 S, and MnS. If a relatively large proportion of 
condensed particles is retained in the atmosphere, cloud layers could affect the structure of the brown dwarf, 
making it hotter by as much as 100 K at 1 bar (depending on the uncertain particle sizes and optical 
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properties). Clouds might increase the flux in the K band, due to the higher temperatures, and lower the 
flux below 1.3 |«m, due to scattering. 

Given these uncertainties, our best fits for the bolometric luminosity, the observed spectrum, and the 
photometry give combinations of T c g and g lying in the range 850 < T e g (K) < 1100 and g < 3000 ms -2 (Fig. 
2). Lower T c g are allowed for g < 300 ms -2 , but the shapes of the J and H bands increasingly deviate from 
the observations. The high- Teg- limit arises from the inability to fit simultaneously the bolometric luminosity 
and the 10 /zm flux. 

A determination of Gl 229 B's gravity via spectral matching would impose a direct constraint on its 
mass. Although g is a function of both mass and radius, the radii of brown dwarfs in this temperature 
range vary relatively little as the mass varies by an order of magnitude. However, at the present stage of the 
analysis the gravity is poorly constrained since high g, high T e ff models fit the spectra as well as lower g and 
T ff ones. The model spectra suggest that high spectral resolution (A/AA> 1000) observations at 1.8-2.1 /jm 
may provide a tighter constraint on g. 

The depth at which the atmosphere becomes convective depends upon the specified model gravity and 
effective temperature. At the highest-pressure point of each model atmosphere, where the temperature- 
pressure profile merges with an adiabat, the interior entropy is calculated for the purpose of matching an 
interior temperature distribution to the given values of (T c s,g). The full evolutionary behavior of a brown 
dwarf is obtained by supplementing previous boundary conditions for objects with masses ~ 0.3 — 15 Mj 
(Jupiter mass units) (13, 14) with our grid of nongrcy model atmospheres. Such evolutionary models are 
needed because R varies with mass and age by up to 30%. The precise radius of the object is impor- 
tant because we must match not only Gl 229 B's spectrum, but also the inferred bolometric luminosity: 
L = AnR 2 aT^ s (a is the Stefan-Boltzmann constant). Our results can be summarized by the following 
approximate fitting formulas (g in m s~ 2 , T e ff in K): 



M = 36 Mj (g/1000) 64 (T cff /1000) 



0.23 



(1) 



* = 1.1 Gyr (s/lOOO) 17 (T cff /1000) 



-2.8 



(2) 



R= 67200 km (.g/1000) 



-0.18 



(Teff/1000) 



0.11 



(3) 
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The effective temperature and surface gravity of Gl 229 B can now be constrained by three sets of 
observations (which are not independent of each other): (i) the observed spectrum from 1 to 2.5 /mi; (ii) the 
broadband flux in several bandpasses from 2 to 13 /mi (10); and (hi) the bolometric luminosity of the object 
(3). These constraints then limit g < 2200 ms -2 and T off = 960 ± 70 K (Fig. 2). Since the reported age of 
Gl 229 A is > 1 Gyr (1), g is further constrained to lie in the range 800 to 2200msec~ 2 (Fig. 3). 

In the atmospheres of Gl 229 B and Jupiter, convection commences as the optical depth to thermal 
photons becomes large, and the temperature profile closely approaches an adiabatic profile at deeper levels 
owing to efficient convection (Fig. 4). In some models, particularly the lower gravity models and those with 
T ff < 900 K, the radiative-equilibrium lapse rate exceeds the adiabatic lapse rate over a several-bar region 
near 1 bar. These atmospheres exhibit two convective regions, a lower region, presumably continuing to 
great depth, and an upper, detached convective zone. Such a detached convective zone is also predicted for 
the atmosphere of Jupiter (15). 

A stellar evolution code and atmosphere models have allowed us to estimate the physical properties of 
the brown dwarf, Gl 229 B. We derive an effective temperature of 960 ± 70 K and a gravity between 800 
and 2200 m s~ 2 . These results translate into masses and ages of 30-55 Mj and 1-5 Gyr, respectively. As 
Eq 1 and Fig. 3 indicate, gravity maps almost directly into mass, and ambiguity in the former results in 
uncertainty in the latter. Since the inferred mass of Gl 229 B exceeds that required for deuterium burning 
(14), deuterium-bearing molecules should not be present in its atmosphere. While the near infrared spectrum 
of Gl 229 B is dominated by H2O, we confirm the presence of CH4 in the atmosphere from our modeling 
of its features at 1.6-1.8 /mi, 2.2-2.4 /mi, and 3.2-3.6 /mi. In addition, we find a flux enhancement in the 
window at 4-5 /mi throughout the T off range from 124 K (Jupiter) through 1300 K, and, hence, that this 
band is a universal diagnostic for brown dwarfs and planets. 



5 



References and Notes 

1. Nakajima, T. et al. Nature 378, 463 (1995). 

2. Oppenheimcr, B.R., Kulkarni, S.R., Matthews, K., & Nakajima, T. Science 270, 1478 (1995). 

3. Matthews, K., Nakajima, T., Kulkarni, S.R. and Oppenheimer B.R. Astrophys. J., submitted. Geballe, 
T.R., Kulkarni, S.R., Woodward, C.E., and Sloan, G.C. Astrophys. J. Lett., in press. 

4. Burrows, A., Hubbard, W. B., Saumon, D., & Luninc, J. I. Astrophys. J. 406, 158 (1993); Lunine, J.I., 
Hubbard, W.B. and Marley, M.S. Astrophys. J., 310, 238 (1986). 

5. Mayor, M. & Queloz, D. Nature, 378, 355 (1995); Marcy, G.W., & Butler, R.P. Astrophys. J. Lett., 
submitted (1996); Butler, R.P. & Marcy, G.W. Astrophys. J. Lett., submitted (1996); Gatewood, G. 
BAAS, in press. 

6. McKay, CP, Pollack, J.B., & Courtin, R. Lcarus 80, 23 (1989); Marley, M.S., McKay, CP., & Pollack, 
J.B. Icarus, submitted. 

7. Anders, E., & Grevesse, N. Geochim. Cosmochim. Acta 53, 197 (1989). 

8. The opacity calculations include collision-induced absorption by H2-H2 [Borysow, A. & Frommhold, 
L. Astrophys. J. 348, L41 (1990)] and H 2 -Hc [Zheng, C. & Borysow, A. Icarus 113, 84 (1995)] and 
references therein, free-free absorption by H^ [Bell, K.L. J. Phys. B 13, 1859 (1980)], bound-free 
absorption by H~ [John, T.L. Astron. & Astrophys. 193, 189 (1988)], and Rayleigh scattering. The 
absorptions of NH 3 , CH4, and PH 3 were calculated using the HITRAN data base [Hilico, J.C, Loete, 
M., & Brown, L.R., Jr. J. of Mol. Spectr. 152, 229 (1992)] with corrections and extensions. Additional 
tabulations [Strong, K., Taylor, F.W., Calcutt, S.B., Remedios, J.J., & Ballard, J. J. Quant. Spectr. 
Radiat. Transfer 50, 363 (1993)] were used where necessary for CH4, especially shortwards of 1.6 [im. 
Data for H 2 and H 2 S were computed from a direct numerical diagonalization [Wattson, R.B., and 
Rothman L.S. JQSRT 48, 763 (1992)] by R.B. Wattson (personal communication). Absorption by CO 
[Pollack et al. Icarus 103, 1 (1993)] and PH 3 opacity was included in the spectral models, but not in 
the temperature profile computation. The baseline models assume the atmosphere to be free of clouds. 

9. For the temperature profile computation, molecular opacity was treated using the k-coefficient method 
[Goody, R., West, R., Chen, L., & Crisp, D. J. Quant. Spectr. Radiat. Transfer 42, 539 (1989)]. 

6 



After a radiative-equilibrium temperature profile was found, the atmosphere was iteratively adjusted to 
self-consistently solve for the size of the convection zones, given the specified internal heat flux. Given 
the radiative-convective temperature-pressure profiles, high-resolution synthetic spectra were generated 
by solving the radiative transfer equation [Bergeron, P., Wesemael, F., and Fontaine, G. Ap.J. 367, 253 
(1991)] Eighteen thousand frequency points were used in the 1 - 15.4 ^m spectral region. These spectra 
were smoothed with a Gaussian-bandpass filter giving a final resolution of A/AA = 600. 

10. Matthews, K., Nakajima, T., Kulkarni, S., & Oppenheimer, B. IAU Circ. #6280 (1995). 

11. Fegley, B., Jr & Loddcrs, K. Icarus 110, 117 (1994). 

12. K. Lodders (personal communication). 

13. Burrows, A., Saumon, D., Guillot, T., Hubbard, W.B., & Lunine, J.I. Nature 375, 299 (1995). 

14. Saumon, D., et al. Astrophys. J. 460, 993 (1996). 

15. Guillot, T., Gautier, D., Chabrier, G., k Mosscr, B. Icarus 112, 337 (1994). 

16. Lindal, G. Astron.J. 103, 967 (1992). 

17. D.S. is a Hubble Fellow. T.G. is supported by the European Space Agency. This research was sup- 
ported by grants from the National Aeronautics & Space Administration, and from the National Science 
Foundation. We thank T. Geballe for digital versions of the Gl 229 B spectrum, K. Lodders for chemical- 
equilibrium calculations, and K. Zahnle for an insightful review. 



7 



Figure Captions 

Fig. 1 A): Synthetic spectra for (bottom to top) T off = 890, 960, 1030 K and g = 1000 m s~ 2 , together with 
data (3) (red line). The three curves in B) are calculated for the same values of T e fj and g; colored 
boxes show the photometric measurements with bandpasses indicated by their width. The red region 
shows the la error on the measurements while yellow gives the 2a error. Yellow triangles show upper 
limits to narrow band fluxes. In both panels, spectral intervals are labeled with the molecules primarily 
responsible for the opacity in that interval. 

Fig. 2 Limits on T c g and g of Gl 229 B. The grey shaded area delimits the effective temperature and gravity 
of model objects which match within 2a the observed bolometric luminosity (3) of Gl 229 B at any age. 
The other areas show limits from fitting the 1 — 2.5 /xm spectrum (vertical lines) and the 2.5 — 13 /jm 
photometry (horizontal lines). 

Fig. 3 The grey shaded area shows the region of overlap of the three constraints from Fig. 2 (the cutoff at low 
g is arbitrary). Solid lines depict the evolution of T e ff and g as various mass brown dwarfs cool. Several 
contours of constant radius (long-dashed curves) and constant age (short-dashed curves) are also shown. 

Fig. 4 Calculated atmospheric structure for Gl 229 B; the dashed curve shows an adiabat corresponding to 
the deep interior temperature profile. For comparison, a profile for Jupiter (16) is shown, along with its 
calculated prolongation into the adiabatic deep interior (dashed curve). 
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